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ABSTRACT 

^ , We observed the nearby, low-density globular cluster M71 (NGC 6838) with 
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the Chandra X-ray Observatory to study its faint X-ray populations. Five X-ray 
sources were found inside the cluster core radius, including the known eclipsing 
binary millisecond pulsar (MSP) PSR J1953+1846A. The X-ray light curve of 
the source coincident with this MSP shows marginal evidence for periodicity at 
the binary period of 4.2 h. Its hard X-ray spectrum and luminosity resemble 
those of other eclipsing binary MSPs in 47 Tuc, suggesting a similar shock origin 
of the X-ray emission. A further 24 X-ray sources were found within the half- 

l> 

associated with M71, somewhat more than predicted, and that 11 ± 6 are back- 
OO , 

ground sources, both galactic and extragalactic. M71 appears to have more X-ray 



mass radius, reaching to a limiting luminosity of 1.5 x 10 30 ergs s 1 (0.3-8 keV). 
From a radial distribution analysis, we find that 18 ± 6 of these 29 sources are 
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sources between Lx = 10 30 -10 31 ergs s _1 than expected by extrapolating from 
other studied clusters using either mass or collision frequency. We explore the 
spectra and variability of these sources, and describe the results of ground-based 
optical counterpart searches. 

Subject headings: globular clusters: individual (M71, NGC 6838) - pulsars: indi- 
vidual (PSR J1953+1846A, M71A) - X-rays: stars, binaries 



1. Introduction 

Globular cluster X-ray sources are of interest for many reasons. Dense globular clusters 
bring stars into close dynamical encounters that lead to the production of X-ray binaries 
(e.g. Hut, Murphy & Verbunt 1991), and studies of globular clusters with different struc- 
tural parameters can elucidate the details of these mechanisms. Globular clusters provide 
concentrations of faint X-ray sources for study of X-ray populations at a known distance, age 
and metallicity. This is of interest even in the least dense globular clusters, where dynam- 
ically formed X-ray binaries may be fewer than those descended from primordial binaries. 
Observations of globular clusters may provide science unique to individual X-ray sources that 
is enabled by knowledge of the X-ray source's distance, reddening, and other properties. An 
ensemble of such results for many clusters may shed light on the evolution of globular clusters 
and their binary populations. 

Bright X-ray sources, associated with accreting neutron stars, have long been under- 
stood to be produced in globular clusters (Clark 1975). Fainter X-ray sources, composed of 
combinations of accreting neutron stars in quiescence, cataclysmic variables (CVs), millisec- 
ond radio pulsars (MSPs), and chromospherically active binaries (ABs), were known in the 
1980s (Hertz & Grindlay 1983) and have been resolved with Chandra's high spatial resolu- 
tion (e.g. Grindlay et al. 2001a, Verbunt & Lewin 2006). A number of globular clusters have 
been observed with Chandra to fairly low X-ray luminosities (< 10 32 ergs s" 1 ). Studies of 
dense clusters have identified large numbers of accreting neutron stars, CVs, ABs and MSPs 
(e.g. Grindlay et al. 2001a, 2001b, Pooley et al. 2002, Becker et al. 2003). Two relatively 
sparse nearby clusters have been carefully studied with Chandra: M4 (Bassa et al. 2004) and 
NGC 288 (Kong et al. 2006). Those studies have indicated that ABs are prevalent in sparse 
clusters, and that CVs are few in number, but possibly larger in number than predicted by 
empirical extrapolations with the density and mass of the cluster core. 

M71 (NGC 6838) is of particular interest due to its close proximity to Earth (4 kpc). A 
moderately short Chandra study can identify X-ray sources down to a few 10 30 ergs s _1 (0.5- 
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2.5 keV), probing the populations of faint CVs, MSPs, and ABs. This cluster is of moderate 
central density (p c = lO 3O5 L /pc 3 ) and shows no evidence for core collapse (its central con- 
centration parameter (Djorgovski 1993) is 1.15). It is moderately reddened (E(B-V)=0.25) 
and has a globular cluster metallicity ([Fe/H]=-0.73, Harris 1996, updated 2003) slightly 
higher than average. The core, half-mass radius, and tidal radii are r c = 0.63', r^ = 1.65', 
and r t = 8.96', respectively (Harris 1996, updated 2003). Neither Einstein (Hertz & Grind- 
lay 1983) nor ROSAT (Verbunt 2001) detected X-ray sources obviously associated with the 
cluster, with the 0.5-2.5 keV ROSAT HRI upper limit being 2.6 x 10 31 ergs s _1 (for an 
assumed bremsstrahlung spectrum with temperature 0.9 keV). The ROSAT HRI did detect 
10 sources in its field of view, but none of these were within the cluster half-mass radius 
(Verbunt 2001). 

Ransom et al. (2003, 2005) and Hessels et al. (2007) reported the presence of a binary 
MSP in M71, with Ppsr = 4.89 ms and P or b = 4.24 h, and the presence at 20 cm of eclipses 
which last roughly 20% of the orbital period. 

In this paper, we report on Chandra X-ray Observatory observations of M71. We discuss 
the observations, data processing, and source detection, providing source lists, in §2J and 
the radial distribution of sources with respect to the nominal cluster center in £|3j In $U 
we present the results of power-law spectral fits to the brighter detected X-ray sources, 
and in §5] we present and discuss the X-ray color-color diagram for these same sources. We 
then discuss candidate counterparts in §5J including those extracted from the 2MASS, USNO 
B1.0, and TYCHO-2 catalogs ( §6.11) . a recent tabulation of variable stars ( §6.3p . X-ray source 
catalogs ( §6.4p . and the MSP in M71 ( §6.51) . Comparison with Hubble Space Telescope (HST) 
observations of M71 ( §6.2[) will be reported in Huang et al. (2008, in preparation). In JT]we 
discuss spectral and temporal properties of the brightest detected sources. A summary of 
our results is given in §SJ 



2. Observations, Data Reprocessing, and Source Detection 

We obtained a 52.4-ks Chandra observation (ObsID 5434) of M71 (nominal center of 
cluster at J2000 RA 19 ft 53 m 46.1 s , DEC 18° 46' 42") on 20-21 December 2004 using the 
Advanced CCD Imaging Spectrometer (ACIS) in very faint (VF), timed-exposure mode, 
with 3.141-s frame time. Starting from the standard Chandra X-ray Center (CXC) pro- 
cessing (ASCDS version number 7.6.7.1 and CALDB 3.2.1; third reprocessing) level 1 files, 
we reprocessed the data without applying pixel randomization. The reprocessing included 
applying the current charge transfer inefficiency correction; selection of the standard ASCA 
grades 0, 2, 3, 4 and 6; and application of the good time filter. We column-cleaned the data 
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using a variant of a method developed at Pennsylvania State Universitjo As is appropriate 
for data taken in VF mode, we used CLEAN55g to reduce the background and remove cos- 
mic ray afterglows. In analyzing data and unless otherwise specified, we utilized events in 
pulse-invariant channels corresponding to 0.3 to 8 keV. 

We searched for X-ray sources in the observed field employing techniques described in 
Tennant (2006) which use a circular-Gaussian approximation to the point spread function 
(PSF). That author gives a rather detailed account of the method in an Appendix, and 
there shows that for the Chandra Deep Field-North (Brandt et al. 2001) it gives results 
consistent with those obtained using the CIAO script wavdetect. Using a 30 ks subset 
of the much longer Deep Field-North observation made it possible to know exactly which 
sources should be found. The method has been used previously in other published work 
(e.g. Chandra observations of the globular cluster M28 reported in Becker et al. 2003). As 
a result of operating in VF mode, background levels were low throughout the observation 
(e.g. ~ 1.4 x 10~ 7 counts/s/pixel over 0.3-8.0 keV for tmix < r^). Therefore within twice the 
M71 half-mass radius (rn = 1.65'), we set the signal-to-noise threshold, S/N, for detection 
to 2.0, but also required the number of source counts to be at least 5 times the statistical 
uncertainty in the local background estimate. The empirical relation derived by Tennant, 
Cmin = (S/N) 2 / 0.81 then implies a point-source sensitivity limit of about 4.9 counts for 
?" Am < 2r/j and in the energy band 0.3-8.0 keV. Because of the increase in PSF size with off- 
axis distance and the associated increase in background within a detection cell, for Tmh > 277, 
we set the S/N threshold for detection to 2.4, and again required the number of source counts 
to be at least 5 times the statistical uncertainty in the local background estimate. The 
point-source sensitivity limit rises to about 7 counts. Following Tennant (2006) we expect a 
completeness limit of about 10 counts. We discuss the luminosities corresponding to these 
count sensitivity limits in £|3] and $11 

The positions of X-ray sources found in this manner inside 2r^ are listed in Table (U 
which is divided into those inside rv, and those with rv, < Tmix < 27V The positions of both 
sets of sources are displayed in Figure [H As can be seen in Figure [TJ the circle with radius 
2r/j extends slightly off the S3 CCD toward S2, which is a front-side illuminated CCD. We 
find one source, s55, inside 2r^ on S2 and include it in Table [TJ In this table, columns 2-5 
give right ascension RA (J2000), declination DEC (J2000), detect cell radius, r ext , and the 
approximate number of X-ray source counts, C x , detected from the source in the 0.3-8.0 keV 
band. Column 6 lists the radius about the source position for inclusion of the source with 



1 See http://www.astro.psu.edu/xray/acis/recipes/index.html 



See http://cxc.harvard.edu/cont-soft/software/clean55. 1.0. html 
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68% confidence. The corresponding radii for inclusion with 95% or 99% confidence, r^ and 
rg 9 , are found by multiplying column 6 by 1.62 or 2.01, respectively. Column 7 gives the 
distance, Tm71, of the source from the nominal center of M71. The remaining columns have 
to do with candidate counterparts and are discussed in §6.11 

Uncertainties in the aspect solution for ACIS-qj imply <7 sys ~ 0"2 (radial, ~ 0.4" at 
90% confidence); to be conservative, we set a sys = 0"2 per axis. A re- registration (boresight) 
analysis using the positions of the 18 2MASS candidate counterparts listed in Table[T]leads to 
the following conclusions: (a) the existing aspect solution leads to a statistically acceptable 
fit with x 2 — 40.4 for 36 degrees of freedom (1 for each axis per counterpart); and (b) 
the best-fit changes in pointings position and roll angle are small (Aa = 0.03 ± 0.07", 
AS = 0.03 ± 0.07", and A6> = —2.3 ± 2.0') and, using the f-test, statistically insignificant. 
Including the candidate counterpart for the MSP in M71 in this analysis does not change 
these conclusions. 

The positions of X-ray sources detected outside 2rh are listed in Table [2j In this table, 
column 7 now lists the CCD on which the source was detected. The CCD S3 is back- 
illuminated and therefore its response extends to lower energies than does the response of 
the other CCDs which are front-illuminated. Due to dither, sources detected close to CCD 
boundaries may have counts on two CCDs as indicated in Table [2j The sources ss33 and 
ss34 are separated by ~ 6", which is significantly larger than the corresponding values for 
rg 9 so we conclude that these are two different sources. 



3. Distributions 

M71 lies near the plane of the Galaxy, with galactic longitude 56.74° and latitude 
—4.56°, suggesting a significant contribution to the detected X-ray sources from galactic 
field sources as well as a contribution from extragalactic sources. In order to determine 
the radial distribution of the X-ray sources detected in the S3 field, we follow the standard 
analysis by adopting the generalized King model profile (see, e.g., Lugger, Cohn, & Grindlay 
1995; Grindlay et al. 2002; Heinke et al. 2006; Lugger et al. 2007) for the projected radial 
distribution function, s(r), given by 

■< r > = c « + irrvkr' (1) 



See http://cxc.harvard.edu/cal/ASPEC T/celmon71 upper left panel of the first figure 
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where cq is the number of background field X-ray sources (galactic and extragalactic) per 
unit solid angle on the sky, and s the number of globular cluster X-ray sources per unit solid 
angle at the cluster center. We make the usual assumption that the projected distribution 
for the visible stars is given by Eq. (1) with 13=1 and r equal to the core radius of 0.63'. 
We also assume the globular cluster X-ray sources are in thermal equilibrium with the stars. 
Then the X-ray core radius of X-ray sources with mass M x = qMo, where Mq is the nominal 
mass of the visible stars, is given by 



( 2 i/A, _ l)VV , (2) 



with p x = (3q - l)/2. 



In contrast to the previous analyses cited above, where the background source level 
was assumed to be entirely of extragalactic origin and thus known a priori, we treated the 
background level as a fitting parameter on an equal basis with the other parameters in 
Eq. (1). In order to determine the best- fit values for c 0jX , s 0tX , and /3 X , we carried out a 
maximum-likelihood fit of Eq. (1) to the radial distances of the sources from the center of 
M71 given in Tabled], but excluded the source on S2. We followed the procedure described 
in Grindlay et al. (2002), which directly fits the unbinned radial distribution, using nonlinear 
optimization to maximize the likelihood. We estimated the parameter value errors using the 
bootstrap method, by generating and fitting 1000 random resamplings, with replacement, of 
the source radial distribution. We took the equivalent 1-er error estimate for each parameter 
to be one half of the 68% range about the median of the distribution of its fitted values from 
the bootstrap resamplings. 

We note that the bootstrap method accounts for Poisson errors in the numbers of both 
cluster and background sources, since each resampling represents a particular realization 
of the underlying probability distribution defined by the original source sample. The total 
number of sources included with any region of the cluster varies over the distribution of 
bootstrap resamplings, since an individual source may be included in each resampling 0, 1, 
2, or more times. For a distribution of 1000 resamplings of the 62-source sample within 
2r^, we find that the mean and standard deviation of the number of sources in the core is 
5.0 ± 2.2, in close agreement with the expectation from Poisson statistics. Thus the King- 
model parameter error estimates produced by the bootstrap method include the contribution 
from Poisson error. 

We determined best-fit parameter values using all S3 sources in Table [T] (i.e. all S3 
sources within 2r^ of the cluster center) and also using just the subset of these sources with 
C x > 10. Our choice of 2r^ for the outer radius of the fit, rather than the more common 
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choice of rh, was based on experiments with a range of values for the outer radius. We 
found that the relatively small ratio of rh/r c = 2.6 for M71 required that we adopt an outer 
radius of at least 2r^ in order to stably determine both the surface density slope, /3 X , and 
the background level, c 0>x . In comparison, 47 Tuc, for which Heinke et al. (2005) adopted r h 
as the outer radius of the fit, has a ratio of rh/r c = 7.0. In performing the fits, we corrected 
for the portion of the circular region (r < 2rh) that lies beyond the edge of S3. The best-fit 
parameters are given in Table [3] and the best-fit model for the radial profile is shown in 
Figure [21 This figure shows the fit to the observed cumulative radial profile, together with 
the separate cluster and background components of the model. 

The large parameter uncertainties listed in Table [3] are the consequence of the high back- 
ground source level, relative to the size of the cluster source population. Table @] compares 
the number of X-ray sources actually detected on CCD S3 within the radii r c , rh, and 2r^ 
with the number of background (field) X-ray sources predicted by the best-fit generalized 
King model. The two-sided errors on these numbers were calculated by adding in quadrature 
the errors propagated from the best-fit and small number Poisson errors using Eqs. (7) and 
(11) from Gehrels (1986). We average the two-sided errors when they differ by < 20%. This 
results in substantial uncertainties for all parameters of the cluster source population, e.g., 
the size of the cluster source population within r h is N x = 18 ± 6. 

Examination of Table [3] indicates that for the entire set of X-ray sources the best-fit 
value of the mass ratio is q — 1.08 ± 0.40, suggesting that the masses of the sources are 
typically of order the mass of the visible stars that dominate the potential (i.e., the turnoff- 
mass stars). However, we note that the large uncertainty in this parameter precludes any 
definitive conclusions about the typical source mass. Similarly, the subsample of brighter 
X-ray sources has a larger best-fit value for the mass ratio, q = 1.36 ± 0.47, but the two q 
values do not differ at a statistically significant level, precluding a test of the dependence of 
source mass on luminosity. 

We may estimate the extragalactic contribution to the background sources as follows. 
Given the value for Eb-v from Harris (1996), the relationships Ay = 3.1Eb-v (Rieke & 
Lebofsky 1985) and n H /A v = 1.79 x 10 21 (Predehl & Schmitt 1995) imply a value for the 
hydrogen column density to M71 of n# = 1.39 x 10 21 H atoms cm -2 . Assuming this value for 
Tin and a power-law spectrum with photon index 1.53 from the Chandra Deep Field-South 
(CDF-S; Giacconi et al. 2001), our detection limit of 4.9 counts in the 0.3-8.0 keV band 
corresponds to a 0.3-8.0 keV flux of 1.32 x 10~ 15 ergs cm -2 s _1 and a luminosity at the 
distance of M71 of 2.5 x 10 30 ergs s -1 . The corresponding value for the 0.5-2.0 keV flux is 
4.2 x 10~ 16 ergs cm" 2 s -1 . Then from the CDF-S (Giacconi et al. 2001, Eq. (1)), we find 
for the expectation value for the number of extragalactic sources above a 0.3-8.0 keV band 



detection limit of 4.9 counts of (0.39 ± 0.09) r\, ni arcmin -2 inside a radius Tmix (in arcmin). 
Repeating this exercise for a detection limit of 10 counts in the 0.3-8.0 keV band leads to 
an expectation value for the number of extragalactic sources of (0.21 ± 0.05) r 2 /71 arcmin -2 
inside a radius rMn (in arcmin). Assuming the same value for nn but using the power-law 
index 1.4 from the Chandra Deep Field-North (CDF-N; Brandt et al. 2001) increases these 
detection limits by 3%. 

Tabled] also lists the number of extragalactic X-ray sources inside the radii r c , r h , and 
2r/j predicted by these results. The errors quoted in the table were again calculated by 
adding in quadrature the expectation value errors and small number Poisson errors using 
Eqs. (7) and (11) from Gehrels (1986), resulting as before in substantial uncertainties. 
Because of the large uncertainties, we cannot conclusively say whether Galactic or extra- 
galactic X-ray sources dominate the sources not associated with M71. Probably both are 
important. Figure [3] shows the distribution of the number of sources with 0.3-8.0 keV band 
counts greater than C, N(> C) vs C for the sources detected inside 2r h . For C < 50 
the best unweighted-least-squares fit power-law index for this distribution is -0.68, while for 
50 < C < 110 the corresponding value for this index is -2.64. 



4. Power-law spectral fits 

In order to estimate X-ray luminosities for the brighter sources within 2rh, we carried 
out fits in XSPEC (Arnaud 1996) to power-law spectra for sources with at least 45 source 
counts in the energy band 0.3-8.0 keV plus the candidate counterpart for the MSP in M71 
with 37.5 source counts, with the results shown in Table [5J For each source we constructed 
response files using the CIAO tool mkacisrmf, and we extracted source and local background 
spectral files using lextrct (Tennant 2006). In all cases we fixed the hydrogen column density 
at the value, nn = 1-39 x 10 21 cm -2 , appropriate for M71. In this table, column 1 gives 
the Chandra source name from Table [1] (see also Figure [1]), column 2 lists the minimum 
number of counts used to group the spectral data for fitting in XSPEC, and columns 3-5 
give the best-fit power-law index, the best-fit normalization, and the value obtained for x 2 
together with the number of degrees of freedom, v. Whenever x 2 l v < 2, we provide the 
single parameter 67% confidence errors calculated by XSPEC using the error command with 
Ax 2 = 1; when x 2 l u > 2, as is the case for source s52, we do not quote any errors. Letting 
the hydrogen column density vary for s52 produces an improved but still unacceptable fit. 
Note that the fits for the sources ssOS ands ss45, while satisfying x 2 l v < 2, are not acceptable 
with better than 99.9% confidence. Column 6 lists the probability of obtaining by chance 
that value for x 2 or greater. Column 7 gives the corresponding unabsorbed X-ray luminosity 
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in the bands 0.5-2.5 keV and 0.3-8.0 keV. 

In order to determine X-ray flux and luminosity limits for the fainter sources within r^, 
we divided the sources with C x < 45, but not including the MSP candidate counterpart, into 
two groups, one with counts in the range 15 < C x < 45 (Group 1 with 7 sources), one with 
C x < 15 (Group 2 with 15 sources). We summed the PI spectra for each group, averaged the 
corresponding responses, and fit the results to power-law spectra as above, with the results 
also shown in Table [5J The total number of counts for the sources in Group 2 (the faintest 
group) is 119.1. Our source count sensitivity limit of 4.9 counts then corresponds to a 0.3-8.0 
keV unabsorbed X-ray flux of ~ 8.0 x 10~ 16 ergs cm -2 s _1 , and a 0.3-8.0 keV unabsorbed 
X-ray luminosity of ~ 1.5 x 10 30 ergs s _1 at the distance of M71. Corresponding limits in the 
0.5-2.5 keV and 0.5-6.0 keV energy bands are ~ 7.9 x 10 29 ergs s _1 and ~ 1.1 x 10 30 ergs s _1 , 
respectively. Very similar limits were found for the somewhat brighter, but fewer, sources in 
Group 1. 

We divided the faint sources with cluster radii between r^ and 2rh into two groups as 
above, Group 3 with 15 < C x < 45 (5 sources) and Group 4 with C x < 15 (22 sources). 
The summed spectra for Group 4 do not fit a power-law spectrum with n# = 1.39 x 10 21 
cm~ 2 particularly well, so we allowed ng to also vary. We then found an acceptable fit to a 
power-law spectrum, with a single parameter 67% upper limit for uh < 2 x 10 20 cm -2 . 



5. X-ray Color— Color Diagram 

Figure H] shows an X-ray color-color diagram for sources inside 2r^, using the three 
energy bands S = C(0.3-0.8 keV), M = C(0.8-2.0 keV), and H = C(2.0-8.0 keV), plotted 
using (H — S)/T as the x-axis and M/T as the y-axis, where T = C(0. 3-8.0 keV). We 
included results for s52 and Group 4 on this plot, even though power-law models with n# 
fixed at the value appropriate for M71 did not provide acceptable fits to the X-ray spectra 
from these sources. Also shown are three curves representing power-law spectra with indices 
ranging from -1 (right end point of curves) to -6 (left end point of curves). The three curves 
are for Hydrogen column densities nn = 0.1, 1.39 (the value for M71), and 10.0 in units of 
10 21 cm" 2 . 

As described by Tennant (2006), using these axes all sources should lie inside the triangle 
defined by S = 0, M = 0, and H = 0. Soft sources will lie to the left, hard sources to the 
right, and centrally peaked sources in the middle. Sources to the left are more likely to be 
stars in our Galaxy, and sources to the right more likely to be pulsars or background AGN. 
Indeed, as we describe in §6.11 s37 and s52, which both have candidate counterparts from 
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the 2MASS and USNO catalogs and are likely stars, lie near the line H = 0, while the MSP 
candidate counterpart s08 lies near the line 5 = 0. The candidate counterpart for s52 is also 
listed in the Tycho-2 catalog and is undoubtedly a foreground star. 



6. Searches for counterparts 

6.1. 2MASS, USNO, & TYCHO-2 catalogs 

Using the HEASARC BROWSE tool in batch mode, we searched the 2MASS9, USNO 
B1.0 (Monet et al. 2003), and TYCHO-2 (Hog et al. 2000) catalogs for possible optical and 
infrared counterparts. The TYCHO-2 catalog is a subset of the USNO B1.0 catalog con- 
taining the 2.5 million brightest stars. Thus coincidence with a TYCHO-2 catalog member 
potentially restricts the nature of the X-ray source. We required that candidates lie inside 
the radius rgg, given in Tables [1] and [2] from the Chandra X-ray source position. The results 
are listed in Tables [TJ for sources with rui\ < 2r^ and Table [2] for sources with 2r^ < tmti, m - 
columns 9 — 11. As noted in §2J we carried out a re- registration (boresight) analysis using the 
X-ray positions and 2MASS postions for the 18 sources in Table [TJ with 2MASS candidate 
counterparts. This analysis demonstrated that for our observation there is no need to shift 
the on-axis pointing position or spacecraft roll angle, thus eliminating the need to search the 
catalogs a second time. Column 8 lists the radial offset from the Chandra position in ". If 
there is a candidate from the 2MASS catalog, column 9 lists the quoted J (the catalog also 
lists H and K magnitudes when available). If there is a candidate from the USNO catalog, 
column 10 lists the Rl magnitude or average of Rl and R2 magnitudes (the catalog also lists 
B and I magnitudes when available). If there is a candidate from the TYCHO-2 catalog, 
column 11 lists the TYCHO-2 VT magnitude. Column 12 lists the catalog and name (if 
there is one) for the candidate counterpart, and in a few cases other pertinent information. 
Column 13 lists the probability of a chance coincidence in %, calculated from 

Bcoinc = JV(> Wl jftrgg / Agearch, (o) 

where N(> m) is the number of sources in the corresponding catalog inside the search area, 
A se arch, with magnitude greater than that of the candidate counterpart, and r 99 the radius 
enclosing the X-ray source location with 99% confidence. For the 2MASS catalog we set 
m = K, and for the USNO B1.0 catalog m — I. The number of possibilities from the 



See http://www.ipac.caltech.edu/2mass/ 
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TYCHO-2 catalog is very small, so P CO inc << 1 in the few cases where a potential TYCHO-2 
candidate counterpart exists. We set the search area, A search , to the area inside the S-array 
boundaries for X-ray sources detected on CCDs S2, S3, or S4, and to the area inside the 
I- array boundaries for X-ray sources detected on CCDs 12 and 13. 

As one might expect, we found no candidates from these catalogs for X-ray sources 
in the core of M71. We found 18 X-ray sources with r c < tmii < 2r/j having a single 
candidate counterpart in the 2MASS catalog, none with multiple counterparts. Of these 18 
X-ray sources, 3 also have a single USNO B1.0 candidate counterpart, and one (s45) has two 
USNO B1.0 candidate counterparts. We found just one X-ray source, s59, having a USNO 
B1.0 candidate counterpart but not a candidate counterpart in the 2MASS catalog. The X- 
ray source s52 has a single counterpart in each of the 2MASS, USNO B1.0, and TYCHO-2 
and is likely a bright foreground star. 

For X-ray sources with 2r^ < tmh, we found 24 had single candidate counterparts in the 
2MASS catalog (20 of these on CCDs S2-S3-S4, 4 on 12-13), none with multiple counterparts. 
Of these 24, 16 also have a single USNO B1.0 candidate counterpart, and three have two 
USNO B1.0 candidate counterparts. There are three X-ray sources with a single USNO B1.0 
candidate counterpart but no 2MASS candidate counterpart, and two X-ray sources with 
two USNO B1.0 candidate counterparts but no 2MASS candidate counterpart. The X-ray 
source ssOS has a single counterpart in each of the 2MASS, USNO B1.0, and TYCHO-2 and 
is likely a bright foreground star. 

Color-color diagrams of J — H vs. J — K, B — RvsB, and R — I vs. R show nothing 
unusual about the candidate counterparts to the Chandra X-ray sources, except that the 
counterparts to s52 and ssOS are both bright in B and R, as might be expected since both 
appear in the TYCHO-2 catalog. 



6.2. HST data 

There are two sets of HST observations of the M71 field currently in the public domain, 
(see HST proposal 8118, Piotto, G., also Piotto et al. 2002; and HST proposal 10524, Ferraro, 
F.). These data consist of 2183 s of WFPC2 data divided among four filters, and 628 s of 
ACS data. Results from comparing these data to the Chandra list of sources presented here 
will be reported in Huang et al. (2008, in preparation). 
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6.3. Variable sources 



We also compared our source positions to those listed for faint variable sources in Table 
9 of Park & Nemec (2000). Of the 23 variable sources listed, 7 lie inside a radius 2r^ of 
the nominal center of M71 but none of them are positionally coincident with a Chandra 
X-ray source position. Two sources from Park & Nemec are coincident with Chandra X-ray 
sources outside 2r/ l , namely ssl4 on CCD S4 and ss52 on CCD S2. The X-ray source ssl4 
lies 0.88" from Park & Nemec source v21, a W UMa binary with period 0.353 d which lies 
in the subgiant region in the Color-Magnitude Diagram. The X-ray source ss52 lies 1.09" 
from Park & Nemec source vl9, which may be an eclipsing binary of thus far unknown type. 
They observed a rise of 0.5 magnitude over ~5 h on 1996 July 12-13 (see Park & Nemec 
2000, Figure 20). Both of these sources have 2MASS and USNO candidate counterparts (see 
Table HJ. 



6.4. X-ray catalogs 

As previously mentioned, Verbunt (2001) reanalyzed all the ROSAT data containing 
M71 and found no X-ray sources inside the cluster's half-mass radius. Again using the 
HEASARC BROWSE tool in batch mode, we searched the HEASARC Master X-ray Catalog 
for sources within 30" of a Chandra X-ray source, with the results listed in Table [61 Inside 
2r h , one source, 1BMW 195344.3+184610, from the Brera Multi-scale Wavelength ROSAT 
High Resolution Imager Catalog (BMWHRICAT, Panzera et al. 2003), lies with within 1" of 
the relatively bright Chandra X-ray source s05 lying just outside the M71 core radius, and 
within 20" or less of the sources s04 and s06. We found no other coincidences inside 2r^. 

Outside 2rh, the bright Chandra X-ray ssOS on CCD S4 is positionally coincident within 
1" of the ROSAT HRI X-ray source 1RXH J195303.5+18520. This source is probably the 
same as other coincidences listed in the table for ss03 drawn from the ROSAT catalogs 
ROSPSPC, WGACAT, and BMWHRICAT. As noted previously, this source is probably a 
bright foreground star. Other coincidences within 30" are listed in the table for the X-ray 
sources ss08, ss45, is06, and, with low probability, isl4- 



6.5. The milli-second pulsar 

The position of PSR J1953+1846A = M71A, the 4.89 ms MSP in a 4.24 h eclipsing 
binary in the core of M71 reported by Ransom et al. (2003, 2005; also see Hessels et al. 2007), 
is separated by 0.12" from the position of the X-ray source s08 in Table [U, well inside the 
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rgg value of 0.65" found for this source. The radio timing data yield a minimum companion 
mass of 0.032 M Q (Ransom et al. 2005, Hessels et al. 2007). As given in Tabled its X-ray 
spectrum is consistent with a power-law spectrum with photon index 7 = 1.89 ± 0.32, and 
an X-ray luminosity at the distance of M71 of 6 x 10 30 ergs s _1 in the 0.5-2.5 keV band 
and 12 x 10 30 ergs s _1 in the 0.3-8.0 keV band. The neutron star atmosphere model (nsa 
in XSPEC; Zavlin et al. 1996, Pavlov et al. 1995) for magnetic field strengths of 0, 10 12 , 
and 10 13 Gauss, with the distance held fixed at 4 kpc, produces bad fits, with the best null 
hypothesis probability being 6.2 x 10 -6 . We therefore conclude that the X-ray emission 
from this MSP is largely nonthermal, and is due either to magnetospheric radiation or to an 
intrabinary shock formed due to interaction between the relativistic pulsar wind and matter 
from the stellar companion (as in the eclipsing MSPs in 47 Tuc, Bogdanov et al. 2006). 

The X-ray luminosity of M71A falls inside the range of X-ray luminosities of the MSPs 
in the well studied globular cluster 47 Tuc (Bogdanov et al. 2006, see their Table 4). It 
also shares the non-thermal nature of its X-ray spectrum with MSPs J, O and W in 47 
Tuc, with a photon power-law index of 7 = 1.89 ± 0.32 vs. ~1 — 1.5 for those three MSPs. 
All four of these binary MSPs show radio eclipses, and PSR J1953+1846A is also eclipsed 
for ~20% of its orbit (Hessels et al. 2007). Bogdanov, Grindlay, & van den Berg (2005) 
reported variability in the X-ray emission from MSP W in 47 Tuc. Setting zero binary phase 
using the radio ephemeris (Stairs et al., in preparation), and folding with 5 phase bins at 
the binary period of 0.1768 d = 4.2431 h = 15.2751 ks leads to the light curve shown in 
Figure O Testing this light curve using Pearson's \ 2 leads to a single trial confidence level of 
2.2% for acceptance of a steady source model, providing marginal (just under 3a) evidence 
in the X-ray band for periodicity at the radio binary period. The phase spanned by the radio 
eclipse (~0.18 — 0.36, also shown in the Figure) does not quite line up with the minimum in 
the X-ray light curve. The error in eclipse phase due to extrapolation from the time of the 
radio observations to the time of the X-ray observation is negligible. For 47 Tuc W Bogdanov 
et al. (2006) found that X-ray minimum precedes the optically derived superior conjunction 
of the pulsar, which may therefore also explain the displacement of the radio eclipse from the 
minimum in the folded X-ray light curve. Bogdanov et al. (2006) suggest the X-ray emission 
for the three eclipsing binary MSPs in 47 Tuc originates in a shock forming from the collision 
of a relativistic wind from the pulsar with material from its low-mass companion, a model 
which is as likely to apply to M71A as to the 47 Tuc sources. 
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Individual sources 



Figure [6] displays source and local background light curves in the 0.3-8.0 keV band for 
the 14 X-ray Chandra detected sources in Tables [1] and [2] with 0.3-8.0 keV counts C x > 80, 
plus the source, s08, in M71's core coincident with the MSP. Four of these sources, s02, s05, 
s08, and s20, lie inside the half- mass radius rv At the timescales shown, most of these 15 
sources are inconsistent with being steady at 99.9% confidence or greater, except for s08 
(the source coincident with the MSP), s52, ss21, and isOl. Source s02 shows a steady rise 
over nearly the last half of the observation. The source ss06 is extremely faint for over half 
of our observation, suddenly flares up strongly, then decays back to a faint level, with the 
flare lasting lasting ~10 ks. This source has a 2MASS candidate counterpart. The strong 
temporal variations in s02 and ss06 suggest flaring on coronally active stars. 

X-ray source s52is positionally coincident with single 2MASS, USNO B1.0, and TYCHO- 
2 sources (see §6.ip . and its Tycho-2 counterpart has significant proper motion. As shown 
in Table its X-ray spectrum fits acceptably to a single MEKAL model with no interstel- 
lar absorption. These facts suggest that s52 is likely a nearby foreground star. Source s38 
shows evidence for time variability (Figure [6]). The USNO catalog provides morphological 
star/galaxy classifications, which allow us to suggest that sources s^5and s59 are stars while 
s37 is a galaxy. 

The brightest of the 73 X-ray sources found outside 2r/ l , indeed the brightest in our field, 
is ss03. This X-ray source is positionally coincident with single 2MASS, USNO B1.0, and 
TYCHO-2 sources (see §6.ip . and on that basis is likely a nearby foreground star. It is also 
coincident with what is likely a single ROSAT X-ray source (see §6.41) . The fit of its X-ray 
spectrum to a power-law model is statistically not very good, but a fit to a single MEKAL 
model is statistically far worse. A fit to a MEKAL+MEKAL model with no interstellar 
absorption is an improvement over a single MEKAL model, but is still statistically worse 
than the fit to a power-law model (compare Tables [5] and [7j). In addition, the source is 
time variable, showing a definite decline during our observation. A single MEKAL model 
does not fit the X-ray spectrum for source ss45, but a MEKAL+MEKAL model with no 
interstellar absorption does fit acceptably. There is one 2MASS object and one USNO B1.0 
object within this source's rgg radius (see §6. 1|) as well as two ROSAT HRI sources within 
2-6". Table [7] shows that the X-ray spectrum for source is06 is fit acceptably by a single 
MEKAL model with no interstellar absorption. We find no counterpart in the 2MASS or 
USNO B1.0 catalogs, but this source is within 2.1" of a faint ROSAT HRI source (see §6.4p . 
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8. Discussion 

We have found five X-ray sources (including the millisecond pulsar M71 A) located within 
the cluster core radius. Our radial distribution analysis indicates 1-2, and possibly all four 
other than M71A, are background sources. M71A shows marginal evidence for modulation 
at the the binary period, slightly before the radio eclipse phase, suggesting similarity with 
the eclipsing millisecond radio pulsar 47 Tuc W (Bogdanov et al. 2006). We identify 29 X- 
ray sources within the half-mass radius of M71, down to a limiting X-ray luminosity (0.3-8.0 
keV) ~ 1.5 x 10 30 ergs s _1 , of which 5-17 may be foreground or background sources. Seven of 
these (all outside the core) have candidate 2MASS counterparts. We also identify a further 
108 sources outside this radius, of which the majority (71) have candidate 2MASS, USNO, 
or TYCHO-2 optical or near-IR counterparts. 

It is of great interest to compare the populations of faint X-ray sources in different 
globular clusters. Our radial distribution analysis shows that the X-ray sources associated 
with M71 seem to be largely confined within the half-mass radius, as seen in denser clusters 
(e.g. NGC 6440, Pooley et al. 2002). This is in contrast to the similarly low-density cluster 
uj Cen, where a quiescent low-mass X-ray binary lies outside the half-mass radius (Rutledge 
et al. 2003). Three likely RS CVn and eclipsing Algol stars (all cluster members) have also 
been astrometrically identified with Chandra or XMM- Newton sources outside the half-mass 
radius (Cool et al. 2002; Gendre, Barret, & Webb 2003). We suggest that this contrast 
may be due to the particularly low level of mass segregation in u Cen, and its much longer 
half-mass relaxation time (in agreement with Verbunt & Johnston 2000). 

Studies of globular clusters show that above Lx = 10 31 ergs s _1 , the X-ray sources seem 
to be dominated by CVs, while below Lx = 10 31 ergs s" 1 , ABs seem to be the dominant 
population (see Fig. 14 in Heinke et al. 2005, or Fig. 6 in Kong et al. 2006). (Neutron stars 
contribute to both populations-as qLMXBs for the brighter population, and as MSPs for 
the fainter-but do not approach 50% of either population in any studied cluster.) Therefore, 
some constraints upon the relative frequency of ABs in different cluster environments are 
possible even before completion of detailed counterpart identification studies. 

We count the total numbers of sources with: a) 10 31 < Lx] and b) 10 30 < Lx < 10 31 
ergs s _1 , in the 0.5-2.5 keV band, inside the half-mass radii of four clusters: M4, M71, 47 
Tucanae, and NGC 6397. We calculate the numbers of background sources within each 
half-mass radius (except NGC 6397, where we use the 2' radius studied by Grindlay et al. 
2001b) predicted by the Chandra Deep Field-North studies (Brandt et al. 2001). For M71, 
we use the counts from Table 1 and derive (0.5-2.5 keV) luminosities using the spectral fit 
to groups 1 and 2 in Table 5, giving Lx = 1-48 x 10 29 ergs/s/photon for on-axis sources. 
M71's high foreground source contamination (attributable to its low Galactic latitude and 
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low density) requires correction for foreground contamination using the results of $3j We find 
a total of 12.5 ± 6.2 X-ray sources within M71's half-mass radius between 10 30 < Lx < 10 31 
ergs/s. The background-subtracted source numbers in these L x ranges are listed in Table El 
and plotted in Figure [7J The errors on the source numbers are derived from the Poisson 
statistics of the numbers of background sources, as in Table HI 

Figure [TJ shows that a linear relation between scaled mass, Mh, inside the half-mass 
radius and the number, N^, of X-ray sources inside the half-mass radius in the luminosity 
range dominated by ABs appears acceptable for most clusters. In computing this linear fit 
(Nh = aMh), we add in quadrature the errors on the predicted number of sources at each 
cluster's mass, with the actual errors on the cluster sources. M71 stands out as having an 
excess of likely ABs for its mass (12.5±6.2) sources with 10 30 < L x < 10 31 ergs s _1 , compared 
to 3.5 expected), although it is less than a 2<r excess. Similar radial profile analyses have 
not yet been done for most of the other comparison clusters; if they also suffer foreground 
contamination, then the discrepancy will be increased. 

The central density of M71 is relatively low, log p c ~ 3.0, suggesting that M71 may 
have lost fewer binaries (and hence ABs) through binary destruction mechanisms (Ivanova 
et al. 2005). Contamination of the sample by collisionally produced X-ray sources (e.g. CVs) 
cannot explain M71's likely abundance of X-ray sources, since M71 is a very low-density and 
low-Gamma cluster (Table 8). Thorough optical identification campaigns (to reduce the 
uncertainties in cluster membership), and deep observations of additional globular clusters, 
will enable progress in understanding the formation and destruction of faint X-ray sources. 
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Table 1. CXO X-ray sources detected with rui\ < 2r^. 
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38.231 
38.252 
38.331 

38.933 
39.648 
39.943 
40.290 
41.397 
41.706 
42.119 
42.438 



19 53 42.835 
19 53 42.953 
19 53 43.794 
19 53 44.536 
19 53 44.967 
19 53 45.904 
19 53 47.022 



19 53 47.711 
19 53 50.931 
19 53 51.292 
19 53 52.201 



18 45 58.21 
18 46 36.74 
18 46 31.92 
18 45 1.67 
18 44 54.57 
18 44 19.88 
18 45 40.38 
18 49 25.52 

18 45 46.61 
18 43 52.84 
18 44 25.46 
18 47 47.23 
18 44 44.55 
18 49 14.34 
18 45 11.14 
18 49 5.27 



18 49 1.01 
18 49 13.41 
18 44 26.65 
18 48 20.00 
18 49 52.03 
18 48 45.91 
18 44 26.78 



18 49 38.07 
18 48 28.42 
18 43 47.55 
18 47 58.22 



2.8 
2.8 
2.4 
2.5 
2.1 
2.3 
1.8 
3.0 

1.7 
2.3 
2.0 
1.8 
1.7 
2.5 
1.5 
2.3 



2.2 
2.4 
1.7 
1.7 
2.8 
1.9 
1.5 



2.5 
1.8 

2.0 
1.6 



fi.O 
6.9 
5.5 
8.6 

11.1 
5.3 
7.1 

45.3 

102.7 

71.1 

04. (i 

32.3 

18.9 

5.8 

4.8 

5.8 



14.8 
6.3 

15.2 
8.6 

31. (i 
5.7 

85.9 



63.4 
27.2 
23.0 

5.0 



0.75 
0.70 
0.68 
0.60 
0.48 
0.66 
0.51 
0.41 

0.32 
0.35 
0.34 
0.36 
0.38 
0.69 
0.51 
0.65 



0.46 
0.65 
0.40 

0.47 
0.43 
0.57 
0.32 



0.36 
0.37 
0.39 

0.53 



3.27 
3.17 
2.71 
3.14 
2.G4 
3.02 
2.13 
3.29 

1.94 
3.21 
2.71 
1.75 
2.26 
2.74 
1.79 
2.54 



2.44 
2.63 
2.33 
1.67 
3.17 
2.06 
2.27 



2.95 
2.11 
3.16 

1.92 



0.18 
0.38 



0.18 
0.98 

0.89 
0.51 

1.10 



0.72 
0.35 



0.01 
0.08 

O.OG 

0.65 



14.8(14.4) 



12.676 
15.751 



16.071 
14.940 



15.1 
14.8(14.5) 



10.5(10.5) 



0.18 
0.32 



2MASS 

USNO 1088-0475990 



0.11 2MASS 
1.51 2MASS 

1.35 2MASS 

USNO 1088-0476353 

1.14 USNO 1088-0476355 



1.28 2MASS 
0.64 2MASS 



0.011 2MASS 
0.035 USNO 1087-0482780 
TYC 1620-1232-1 



0.16 2MASS 
2.38 2MASS 
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Table 1 — Continued 



(i) 



(2) 



(3) 



(4) 



(■ r >) 



(6) 



(7) 



(8) 



(9) 



(1(1) 



(11) (12) 



(13) 



s57 


19 53 52.738 


18 47 52.18 


1.6 


8.8 


0.44 


1.95 


s58 


19 53 53.482 


18 45 29.41 


1.4 


5.3 


0.48 


2.13 


s59 


19 53 54.710 


18 48 38.75 


2.2 


12.7 


0.47 


2.81 


s60 


19 53 56.921 


18 48 16.53 


2.2 


8.5 


0.55 


3.00 


s61 


19 53 58.155 


18 46 49.10 


1.9 


11.0 


0.46 


2.85 


s62 


19 53 58.396 


18 47 18.34 


2.0 


6.1 


0.58 


2.97 


s63 


19 53 59.410 


18 47 26.64 


2.2 


5.2 


0.67 


3.23 



17.8(17.5) 



USNO 1088-0477382 



a R.adius of the detect cell for collecting X-ray counts. 

Number of 0.3—8.0 kcV source counts collected in the detect cell. 

c Radius enclosing the true source position with 68% confidence. The corresponding radii for inclusion with 95% or 99% confidence are 
found by multiplying this column by 1.62 or 2.01, respectively. 

d Radius from the nominal center of the cluster (J2000 RA 19 h 53 m 46.1 s , DEC +18° 46 ; 42"). 

°R.adial offset of candidate counterparts found by searching the HEASAR.C 2MASS (B/2mass; sec http://www.ipac.caltcch.edu/2mass/ I, 
USNO-B1.0 (1/284; Monet et al. 2003), and TYCHO-2 (tycho2; Hog et al. 2000) catalogs, requiring that the candidate counterpart lie within 
the 99% confidence radius of a CXO source. There arc 6,414 and 3,376 2MASS sources, 15,324 and 9,196 USNO sources, and 7 and 9 TYCHO 
sources, appearing within the S2-S3-S4 and 12-13 boundaries shown in Figures 3 and 4, respectively. If blank, no candidate counterpart was 
found in these catalogs. 

.7 magnitude of a candidate counterpart found in the 2MASS catalog (HEASARC B/2mass; see http://www.ipac.caltech.edu/2mass/ I. 

S R1(R2) magnitude of a candidate counterpart found in the USNO-B1.0 catalog (HEASARC 1/284; Monet et al. 2003). 

h VT magnitude of a candidate counterpart found in the TYCHO-2 catalog (HEASARC tycho2; Hog et al. 2000). 

"Probability of chance coincidence (sec jj6.1| . 

J This source is on the S2 front-side CCD. All other sources in this table are on the S3 back-side CCD. 



-23- 



Table 2. CXO X-ray sources detected with 2rh < tmti- 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Source 


RA(J2000) 

h m s 


Dec(J2000) 

o / // 


fnxt* 


c x b 


■?"x C 


Chip d 


Offset*" 


mag 


mag 


VT h 
mag 


Paoinc' 
% 


Comment 



ssOl 


19 52 


56.900 


18 


51 8.75 


22.7 


49.9 


1.97 


S4 








ss02 


19 53 


2.976 


18 


49 23.95 


16.7 


45.2 


1.53 


S4 








ss03 


19 53 


3.490 


18 


51 59.52 


19.4 


7980.2 


0.33 


S4 


0.08 

0.14 
0.13 


8.813 


10.7(10.6) 


ss04 


19 53 


4.399 


18 


48 38.23 


15.4 


38.2 


1.53 


S4 








ss05 


19 53 


9.151 


18 


53 50.29 


19.1 


85.0 


1.29 


S4 


1.02 

1.21 


14.360 


16.2(15.8) 


ss06 


19 53 


10.843 


18 


50 18.08 


13.0 


251.6 


0.58 


S4 


0.20 


13.918 




ss07 


19 53 


12.205 


18 


47 59.83 


10.fi 


32.5 


1.16 


S4 








ss08 


19 53 


14.354 


18 


51 57.83 


13.3 


386.1 


0.51 


S4 








ss09 


19 53 


14.636 


18 


50 23.35 


11.2 


25.7 


1.36 


S4 


1.48 
1.92 


14.577 


17.3(16.6) 


sslO 


19 53 


16.567 


18 46 37.74 


8.1 


59.9 


0.70 


S3 


1.08 


15.345 




ssll 


19 53 


21.991 


18 48 39.87 


6.6 


35.9 


0.73 


S3 


1.00 


14.916 




ssl2 


19 53 


24.427 


18 49 53.59 


6.7 


20.6 


0.94 


S3-S4 








ss!3 


19 53 


25.045 


18 52 43.85 


10.1 


14.8 


1.62 


S4 


1.38 


16.431 




" 


















1.62 




18.5(17.0) 


ssl4 


19 53 


25.618 


18 


51 17.83 


7.9 


14.9 


1.27 


S4 


l.Ofi 
2.44 


14.397 


14.5(13.9) 


ssl5 


19 53 


26.671 


18 


50 40.09 


6.8 


19.2 


0.98 


S3-S4 








sslfi 


19 53 


27.949 


18 


47 4.78 


4.0 


12.fi 


0.74 


S3 


1.16 

1.48 


13.743 


14.7(14.4) 


ssl7 


19 53 


28.475 


18 


49 51.63 


5.4 


14.7 


0.90 


S3 








ssl8 


19 53 


28.617 


18 


46 58.16 


3.8 


89.9 


0.39 


S3 








ssl9 


19 53 


29.129 


18 


47 13.87 


3.7 


16.8 


0.63 


S3 


0.36 

0.49 


12.603 


13.9(13.5) 


ss20 


19 53 


29.353 


18 46 17.54 


3.6 


29.8 


0.50 


S3 








ss21 


19 53 


29.427 


18 46 39.79 


3.6 


98.9 


0.37 


S3 


0.47 


15.358 




ss22 


19 53 


29.494 


18 


50 8.47 


5.4 


12.fi 


0.97 


S3 








ss23 


19 53 


29.996 


18 


45 19.29 


3.5 


7.3 


0.84 


S3 








ss24 


19 53 


30.251 


18 


51 57.63 


7.4 


20.1 


1.04 


S3-S4 


1.17 




19.9(19.9) 


ss25 


19 53 


31.144 


18 47 39.25 


3.3 


12.6 


0.64 


S3 








ss26 


19 53 


31.853 


18 


48 33.84 


3.6 


12.8 


0.68 


S3 


0.06 
0.68 


13.051 


14.9(14.6) 


ss27 


19 53 


35.327 


18 


50 44.09 


4.8 


7.9 


1.06 


S3 








ss28 


19 53 


37.543 


18 


43 57.35 


2.6 


8.9 


0.60 


S3 


0.22 


13.577 




ss29 


19 53 40.997 


18 


43 26.37 


2.5 


28.8 


0.41 


S3 


0.19 


15.379 




ss30 


19 53 41.110 


18 


50 42.30 


3.9 


8.6 


0.86 


S3 


1.39 




17.1(16.9) 


" 


















1.22 






ss31 


19 53 43.012 


18 


42 14.78 


3.4 


14.fi 


0.61 


S3 


0.36 


14.522 




" 


















0.47 




15.2(15.0) 


ss32 


19 53 44.491 


18 


51 28.40 


4.6 


28.9 


0.60 


S3 








t^33 


19 53 46.137 


18 


42 52.38 


2.6 


8.3 


0.63 


S2-S3 


0.53 


14.809 




" 


















0.76 




15.3(15.0) 


ss34 


19 53 46.572 


18 


42 55.24 


2.6 


8.2 


0.62 


S2-S3 








**35 


19 53 48.153 


18 


43 15.76 


2.3 


10.2 


0.53 


S2-S3 








**3(i 


19 53 48.183 


18 


50 7.69 


3.0 


14.8 


0.56 


S3 








ss37 


19 53 


51.772 


18 


40 42.93 


5.2 


9.5 


l.Ofi 


S2 








ss38 


19 53 


52.141 


18 


42 10.21 


3.4 


14.0 


0.63 


S2 








ss39 


19 53 


57.803 


18 


43 20.46 


2.9 


11.2 


O.fil 


S2 


0.79 




(17.7) 


ss40 


19 53 


58.657 


18 


48 26.15 


2.5 


38.6 


0.39 


S3 








ss41 


19 54 


0.169 


18 


39 29.30 


8.1 


44.7 


0.79 


S2 








ss42 


19 54 


0.240 


18 


41 27.64 


5.2 


30.2 


0.65 


S2 


0.40 

1.02 


14.693 


16.1(15.5) 


ss43 


19 54 


0.960 


18 


43 29.97 


3.3 


8.2 


0.77 


S2 








ss44 


19 54 


1.376 


18 


44 27.20 


2.9 


22.8 


0.47 


S2 








ss45 


19 54 


2.852 


18 


42 50.70 


4.3 


797.2 


0.32 


S2 


0.39 

0.62 


10.654 


12.3(12.2) 


ss46 


19 54 


5.729 


18 


40 58.72 


7.0 


18.3 


1.04 


S2 








ss47 


19 54 


6.665 


18 


42 5.32 


5.9 


9.7 


1.19 


S2 








ss48 


19 54 


6.637 


18 


42 42.25 


5.3 


41.6 


0.58 


S2 








ss49 


19 54 


7.447 


18 


47 17.84 


4.0 


17.4 


0.65 


S2 








ss50 


19 54 


7.856 


18 


44 13.20 


4.5 


10.8 


0.88 


S2 









0.0035 2MASS 
0.0028 USNO 1088-0473175 
TYC 1624-1644-1 

3.33 2MASS 

6.33 USNO 1088-0473600 

0.50 2MASS 



3.63 2MASS 

9.20 USNO 1088-0474040 

1.82 2MASS 

1.35 2MASS 

14.41 2MASS 

15.91 USNO 1088-0474864 

3.42 2MASS 

2.05 USNO 1088-0474902 

0.81 2MASS 

1.12 USNO 1087-0481208 



0.36 

0.48 



2MASS 

USNO 1087-0481286 



0.76 2MASS 



6.82 


USNO 1088-0475303 


0.40 


2MASS 


0.78 


USNO 1088-0475453 


0.54 


2MASS 


0.88 


2MASS 


5.40 


USNO 1088-0476227 




USNO 1088-0476239 


0.96 


2MASS 


1.13 


USNO 1087-0482479 


1.09 


2MASS 


1.03 


USNO 1087-0482712 



4.50 USNO 1087-0483662 



1.02 2MASS 

1.41 USNO 1086-0482319 



0.015 2MASS 

0.053 USNO 1087-0484047 
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Table 2 — Continued 



(i) 



(2) 



(3) 



M 



(5) 



(6) (7) 



(8) 



(9) 



(10) 



(11) 



(12) 



(13) 



ss51 


1!) 


54 


9.007 


18 46 


36.24 


4.3 


40.6 


0.51 


S2 


0.95 


ss52 


19 


54 


9.562 


18 47 


22.66 


4.7 


16.7 


0.75 


S2 


0.88 
0.70 


ss. r >3 


111 


54 


10.927 


18 39 


50.05 


10.3 


15.1 


1.62 


S2 




ss54 


1!) 


54 


11.625 


18 46 


47.69 


5.2 


19.3 


0.78 


S2 




SriS5 


1!) 


54 


11.738 


18 43 


12.42 


6.4 


35.6 


0.72 


S2 




ss56 


1!) 


54 


12.698 


18 43 


45.91 


6.4 


21.9 


0.88 


S2 




ss57 


19 


54 


13.904 


18 37 45.12 


15.3 


30.6 


1.69 


S2 


2.7(1 


" 




















2.32 


ss58 


19 


54 


16.297 


18 39 


3.75 


13.6 


72.7 


1.01 


S2 


1.15 


ss59 


19 


54 


18.319 


18 43 


39.81 


8.7 


46.6 


0.82 


S2 




isOl 


19 


52 


54.455 


18 40 


11.34 


26.3 


106.3 


1.57 


13 




is02 


19 


53 


17.086 


18 42 


34.97 


9.7 


25.3 


1.20 


13 




is03 


19 


53 


22.003 


18 38 


42.71 


13.2 


26.6 


1.58 


13 


1.08 
0.31 
2.92 


is04 


19 


53 


24.281 


18 33 


12.71 


27.1 


71.4 


1.96 


12 


3.60 
3.09 


is05 


19 


53 


25.238 


18 39 


14.66 


11.2 


42.0 


1.08 


13 




is06 


19 


53 


26.354 


18 39 


55.13 


9.6 


144.5 


0.57 


13 




is07 


19 


53 


26.619 


18 41 


3.48 


7.8 


26.4 


0.97 


13 




is08 


19 


53 


28.215 


18 40 


53.31 


7.6 


20.7 


1.05 


13 


0.80 
0.20 

1.92 


is09 


19 


53 


28.730 


18 40 


7.81 


8.6 


9.3 


1.73 


13 


3.44 
2.83 


islO 


19 


53 


31.432 


18 39 


40.69 


8.6 


12.3 


1.52 


12 




isll 


19 


53 


34.916 


18 38 


28.26 


10.2 


14.2 


1.65 


12 




isl2 


19 


53 


38.072 


18 37 44.80 


11.2 


13.1 


1.89 


12 




isl3 


19 


53 47.875 


18 33 


16.88 


22.8 


37.9 


2.26 


12 


3.39 


" 




















2.97 


" 




















3.42 


isl4 


1!) 


53 


55.822 


18 34 


51.53 


18.3 


65.7 


1.40 


12 





18.1(18.4) 
18.0(17.7) 



3.46 USNO 1087-0484528 

2.72 2MASS 

5.56 USNO 1087-0484568 





8.33 


2MASS 


19.04 


50.39 


USNO 1086-0483424 


19.9(19.7) 




USNO 1086-0483626 





10.62 


2MASS 


18.1(17.7) 


13.33 


USNO 1086-0479228 


19.5 


35.44 


USNO 1086-0479239 




15.03 


2MASS 


16.5(15.0) 


4.83 


USNO 1085-0482335 





3.75 


2MASS 


15.8(15.4) 


2.47 


USNO 1086-0479720 


18.6 


2.27 


USNO 1086-0479736 


18.2 


38.88 


USNO 1086-0479755 


19.1(19.3) 


37.82 


USNO 1086-0479756 





9.89 


2MASS 


(15.2) 




USNO 1085-0484067 


15.8(15.7) 


10.93 


USNO 1085-0484068 



a R.adius of the detect cell for collecting X-ray counts. 

Number of 0.3—8.0 keV source counts collected in the detect cell. 

c R.adius enclosing the true source position with 68% confidence. The corresponding radii for inclusion with 95% or 99% confidence arc found 
by multiplying this column by 1.62 or 2.01, respectively 

CCD on which the X-ray source appears. In a few cases, the source is dithered across the gap between CCDs. 

c R.adial offset of candidate counterparts found by searching the HEASARC 2MASS (B/2mass; sec http://www.ipac.caltcch.cdu/2mass/l, 
USNO-B1.0 (1/284; Monet et al. 2003), and TYCHO-2 (tycho2; Hog et al. 2000) catalogs, requiring that the candidate counterpart lie within 
the 99% confidence radius of a CXO source. There arc 6,414 and 3,376 2MASS sources, 15,324 and 9,196 USNO sources, and 7 and 9 TYCHO 
sources, appearing within the S2-S3-S4 and 12-13 boundaries shown in Figures 3 and 4, respectively. If blank, no candidate counterpart was 
found in these catalogs. 

J magnitude of a candidate counterpart found in the 2MASS catalog (HEASAR.C B/2mass; sec http://www.ipac.caltcch.cdu/2mass/l. 

S R1(R2) magnitude of a candidate counterpart found in the USNO-B1.0 catalog (HEASARC 1/284; Monet et al. 2003). 

h VT magnitude of a candidate counterpart found in the TYCHO-2 catalog (HEASARC tycho2; Hog et al. 2000). 

Probability of chance coincidence (sec £|6. H . 
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Table 3. Best-fit King model parameters. 



(i) 



(2) 



(3) 



King model parameter All 62 S3 X-ray sources 33 S3 sources with C' x > 10 



qx 

rcA') 



1.28 ±0.60 
6.92 ± 3.42 
1.13 ±0.61 
1.08 ± 0.40 
0.58 ±0.26 



0.63 ±0.31 
7.17 ±3.65 
1.54 ±0.70 
1.36 ±0.47 
0.48 ±0.16 



a Units are sources per square-arcminutc. 



Table 4. Predicted No. of S3 sources inside radius r Mn 



(1) 


(2) 


(3) 


(4) (5) 


Radius a 


No. detected 15 


Predicted no. of 
background sources 


Extragalactic contribution' 1 
CDF-S C CDF-N f 


All X-ray sources 



r c = 0.63' 

r h = 1.65' 

2r h = 3.30' 



5 
29 
62 



1 6+ 2 ' 6 
1 - D -1.4 

11.0 ±6.4 
39.4±19.7 g 



-+2-1 
-0.5 
,+3.1 
-1.9 
10.8±4.6 g 



0.5: 
3.3"* 



6+ 21 

q q+3.2 

3.«_ 20 

12.2±4.7 g 



X-ray sources with C x > 10 



r c = 0.63' 

r h = 1.65' 

2r h = 3.30' 



1 
18 
33 



5-4l 4 3 1 
19.4±10.f 



n o+2.0 

1 H+ 2 - 6 
1- S -1.3 

+3.8 



5.8^ 



2.7 



0.3 
2.3 

7.5 



2.0 

0.3 
+2.8 

1.5 
+4.1g 

3.0 



a Radius in '. 

b No. of X-ray sources detected on CCD S3 inside this radius. 

c Predicted no. of background (field) X-ray sources on S3 inside this radius, based 
on the value for crj^ from Table [3] See the 502 text for discussion of the assigned 
errors. 

d Predicted no. of extragalactic X-ray sources on S3 inside this radius, based 
on results from the Chandra Deep Fields North and South. See the ^3\ text for 
discussion of the assigned errors. 

e Giacconi et al. (2001). 

f Brandt et al. (2001). 

SReduced by the amount of area inside 2r^ that falls off S3. 
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Table 5. Power-law spectral fits for selected X-ray sources. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(8) 


Source 


Grouping a 


njj/10 22 


7 


A b 


x 2 M 


P (> x 2)d 


L x 


Comment 








X-ray sources on S3 with 


rjvf7l < 2r h l 








s02 


15 


0.139 


3.09 ±0.22 


3.01 ±0.34 


5.82/ 3 


0.12 


15 (31) 




s05 


20 


0.139 


1.57 ±0.11 


7.83 ± 0.62 


17.41/12 


0.14 


41 (103) 




s08 


12 


0.139 


1.89 ±0.32 


1.12 ±0.23 


0.80/ 1 


0.37 


6(12) 


MSP 


si 5 


15 


0.139 


2.18 ±0.24 


2.04 ±0.30 


2.10/ 2 


0.35 


10 (19) 




sl9 


15 


0.139 


2.35 ±0.29 


1.95 ±0.30 


0.0011/ 1 


0.97 


9(18) 




s20 


15 


0.139 


2.17 ±0.19 


3.03 ±0.37 


4.91/ 3 


0.18 


15 (29) 




s29 


15 


0.139 


1.44 ±0.26 


1.30 ±0.26 


0.020/1 


0.89 


7(19) 




s37 


15 


0.139 


2.56 ±0.33 


1.59 ±0.26 


1.37/ 1 


0.24 


8(14) 




s38 


20 


0.139 


1.40 ±0.18 


2.44 ±0.34 


2.60/ 2 


0.27 


13 (37) 




s39 


15 


0.139 


1.36 ±0.22 


1.75 ±0.30 


0.16/ 2 


0.92 


10 (28) 




s40 


15 


0.139 


0.20 ±0.21 


0.66 ±0.16 


1.40/ 2 


0.50 


6(47) 




s52 


15 


0.139 


3.69 


2.21 


23.98/ 3 


0.000025 


12 (31) 


X 2 /^>2 


s53 


20 


0.139 


2.87 ±0.28 


2.43 ±0.32 


1.37/ 1 


0.24 


12 (23) 




Summed spectra for faint sources on S3 with r^Tl < ^ r h 


Group l g 


20 


0.139 


1.72 ±0.14 


5.43 ±0.52 


6.37/ 7 


0.50 


28 (64) 




Group 2 h 


20 


0.139 


2.25 ±0.20 


3.99 ±0.43 


2.90/ 4 


0.58 


19 (37) 




Groups 1+2 1 


20 


0.139 


1.89 ±0.11 


9.44 ±0.68 


8.34/13 


0.82 


47 (100) 




Group 3 J 


20 


0.139 


1.74 ±0.16 




3.30/ 4 


0.51 


19 (43) 




Group 4 k 


20 


0.139 


2.04 


5.28 


14.00/ 7 


0.051 


26 (52) 


X 2 /^>2 


Group 4 


20 


< 0.02 1 


1.44 ±0.16 


n o 9 +0. 40 
°*°'*— 0.30 


5.86/ 6 


0.44 


18 (49) 




Groups 3+4 m 


20 


0.139 


1.90 ±0.14 


9.07 ±0.74 


16.30/13 


0.23 


45 (96) 










X-ray sources with 2r/ l 


< »*M7l n 








ss03 


20 


0.24 ±0.013 


2.55 ± 0.047 


601 ± 24 


225.55/169 


0.0014 


2900 (5400) 


on S4 


ss06 


20 


lg +0.060 

u,ly -0.076 


o t-n+0.49 
• 3 -° -0.26 


20.54^-85 


9.12/ 9 


0.43 


110 (260) 


on S4 


ss08 


20 


29+ - 11 
u,zy -0.075 


1 84+ ' 19 


20.79±4;g 


15.97/ 16 


0.46 


100 (230) 


on S4 


ss45 


20 


0.23 ± 0.031 


3.11 ±0.15 


58 ±8 


56.35/ 31 


0.0035 


290 (600) 


on S2 


isOl 


20 


17+ - 25 
u - 1( -0.14 


1 ™+ ' 26 


„ + 4.06 
°-°'-2.16 


1.28/ 5 


0.94 


30 (78) 


on 13 


is06 


20 


0.46 


5.43 


31.27 


12.96/ 4 


0.011 


300 (1700) 


on 13; x 2 jv > 2 



a Minimum number of counts per spectral bin for fitting in XSPEC. 

b Power-law spectra normalization in units of 10 -6 . 

c Value of x 2 for best-fit and the number of degrees of freedom, v. 

d Probability of finding a value for x 2 > the value actually found. 

c Unabsorbed X-ray luminosity in units of 10 30 ergs s _1 assuming the distance to M71; energy bands are 0.5—2.5 (0.3—8.0) keV. 

i dN/dE(E) = e- n n cr 'SM A/E~* in units of photons/s-cm 2 -keV with E in keV, with n H fixed at the value 1.39 X 10 21 H atoms 
per cm 2 . The quoted errors are the one-parameter 67% confidence level errors given by the XSPEC command error with A\ 2 = 1. 
Errors are not quoted when x 2 /v > 2, as is the case for source s52 and Group 4 with rijj fixed. 

s Results of fit for summed power-law spectrum for Group 1, which contains the 7 faint sources on S3 not included in the table but 
with C x > 15 and with tm71 < f'h- s04, s09, sll, sl3, sl7, sl8, and s26. From Table IT1 the number of source counts in this group is 
199.0. 

h Results of power-law spectral fit for summed power-law spectrum for Group 2, which contains the 15 faintest sources with 
rum < r h : sOl, s03, s06, s07, slO, sl2, sl4, sl6, s21, s22, s23, s24, s25, s27, and s28. From Table[TJ the number of source counts in 
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this group is 119.1. 

'Results of fit for summed power-law spectrum for all sources in Groups 1 and 2. The number of source counts in these groups is 
318.1. 

J Results of fit for summed power-law spectrum for Group 3, which contains the 5 faint sources on S3 not included in the table but 
with C x > 15 and with r% < vm71 < 2r^: s41, s42, s48, s50, s54. From Table[TJ the number of source counts in this group is 125.1. 

k Results of fit for summed power-law spectrum for Group 4, which contains the 22 faintest sources with r^ < rji/71 < 2r^: s30, 
831, s32, s33, s34, s35, s36, s43, s44, s45, s46, s47, s49, s51, s56, s57, s58, s59, s60, s61, s62, s63. From Table[T] the number of source 
counts in this group is 165.1. 

'Since the fit for Group 4 with n^j fixed is not very good, we also fit these data with Njj free. This produced an acceptable fit with 
a single parameter 67% upper limit for nu of ~ 2 X 10 20 cm -2 . 

m Results of fit for summed power-law spectrum for all sources in Groups 3 and 4. 

n dN/dE(E) = e~ nHlJISM A/E 1 in units of photons/s-cm 2 -keV with E in keV. The quoted errors are the one-parameter 67% 
confidence level errors given by the XSPEC command error with A\ 2 = 1. Errors are not quoted when X 2 / u ^ 2, as is the case for 
source is06. 
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Table 6. X-ray sources within 30" of a CXO source. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Source 


CCD 


Catalog 3, 


Name 


Offset 


RA(J2000) 

h m s 


Dec(J2000) 

o / It 


10 -3 c/s 



s05 


S3 


ROSAT BMWHRICAT 


s()4 


S3 




s06 


S3 




ss03 


S4 


ROSAT ROSHRI 

ROSAT ROSPSPC 
ROSAT BMWHRICAT 

ROSAT WGACAT 


ss08 


S4 


ROSAT BMWHRICAT 


ss45 


S2 


ROSAT ROSHRI 
ROSAT BMWHRICAT 


is06 


f:-i 


ROSAT BMWHRICAT 


isl4 


f2 


ROSAT BMWHRICAT 



1BMW 195344.3+184010 



1RXH J195303. 5+18520 
1RXH J195303. 2+18520 
2RXP J195303. 4+18515 
1BMW 195303.3+185201 
1BMW 195303.8+185201 
1WGA J1953. 0+1852 
1BMW 195314.4+185157 
1BMW 195314.0+185200 
1RXH J195402. 9+18424 
1RXH J195402. 6+18425 
1BMW 195402.8+184249 
1BMW 195402.5+184252 
1BMW 195326.2+183955 
1BMW 195355.5+183516 



0.88 


19 


53 44.33 


18 46 


9.98 


0.46+0.12 


13.52 


19 


53 44.33 


18 46 


9.98 


0.46+0.12 


19.24 


19 


53 44.33 


18 46 


9.98 


0.46+0.12 


0.91 


19 


53 


3.5 


18 52 


0.41 


10.69+0.73 


5.67 


19 


53 


3.19 


18 52 


3.29 


8.72+0.73 


3.35 


19 


53 


3.34 


18 51 


56.99 


32.42+2.31 


3.69 


19 


53 


3.24 


18 52 


0.52 


8.89+0.56 


4.08 


19 


53 


3.77 


18 52 


0.59 


10.65+0.58 


6.26 


19 


53 


3.29 


18 52 


5.09 


33.6+2.60 


0.61 


19 


53 


14.38 


18 51 


57.31 


2.22+0.29 


5.58 


19 


53 


13.99 


18 52 


0.01 


1.46+0.23 


2.58 


19 


54 


2.93 


18 42 48.35 


1.85+0.36 


3.71 


19 


54 


2.59 


18 42 


50.94 


2.13+0.37 


2.09 


19 


54 


2.81 


18 42 48.71 


1.66+0.23 


5.94 


19 


54 


2.45 


18 42 


52.2 


2.23+0.28 


2.10 


19 


53 


26.21 


18 39 


55.4 


0.79+0.17 


24.44 


19 


53 


55.49 


18 35 


15.5 


1.74+0.35 



a Sourccs in this table were extracted from the HEASARC Master X-ray Catalog, requiring that an entry lie within 30 of a CXO 
source. ROSHRI and ROSPSPC arc the standard catalogs of pointed observations with the ROSAT HRI and PSPC instruments. 
WGACAT is another catalog of PSPC pointed observations analyzed using the XIMAGE tool (White et al. 1995). BMWHRICAT 
is the Brcra Multi-scale Wavelet ROSAT High Resolution Imager Source Catalog (BMW-HRI) (Panzcra et al. 2003), derived from 
all ROSAT HRI pointed observations with exposure time longer than 100 seconds available in the R.OSAT public archives, and 
analyzed using a wavelet detection algorithm. 

Measured count rate. 



Table 7. MEKAL and MEKAL+MEKAL spectral fits. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Source 


Mo del a 


Grouping 


Tl 


A 1 /10 6 


T 2 


A 2 /10 6 


x 2 /" c 


P(>x 2 ) d 


-Fx° 



s52 MEKAL 

ss03 MEKAL+MEKAL 

ss45 MEKAL+MEKAL 
is06 MEKAL 



15 


0.41 ± 0.06 


2.17 + 0.29 






0.81/ 3 


0.85 


10(11) 


20 


4.00 ± 0.21 


79018 


0.96 ± 0.06 


49 ± 12 


322/165 


3.37c-12 


1470(2800) 


20 


0.83 ± 0.06 


9.9 ± 1.8 


2.84 + 0.33 


51 + 4 


28.1/30 


0.56 


130(200) 


20 


1.03 ± 0.08 


8.12 ± 0.90 






6.12/ 5 


0.29 


29(34) 



a The XSPEC MEKAL model calculates the emission from a hot diffuse gas at temperature T and with normalization A as given at 
http://hcasarc.gsfc.nasa.gov/docs/xanadu/xspcc/manual/XSmodclMckal.html (Mcwe, Groncnschild, &i van den Oord 1985; Mewe, Lcmcn, 
& van den Oord 1986; Kaastra 1992; Liedahl, Osterheld, &: Goldstein 1995; Arnaud &; Rothenflug 1985; Arnaud &; Raymond 1992). The 
quoted errors are the one-parameter 67% confidence level errors given by the XSPEC command error with Ax — 1- 

Minimum number of counts per spectral bin for fitting in XSPEC. 

c Value of x f° r best-fit and the number of degrees of freedom, v. 

Probability of finding a value for x 7l the value actually found. 
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Table 8. Background-subtracted source numbers for selected globular clusters. 



(1) 




(2) 


(3) 


(4) 


(5) 






M4 


M71 


47Tuc 


NGC 6397 


Distance (kpc) 




1.73 


1.0 


4.5 


2.55 


pa 




1 


0.11 


24.9 


2.52 


M h h 




1 


0.30 


10.2 


0.78 


10 31 < L x (ergs s~ 


L )C 


l d 


2.2 


23 


8 


10 30 < L^ (ergs s- r ) c 


< 10 31 


6.6 


12.5 


116 


7 


Error 




3.0 


6.2 


3.9 


2.1 


References 




1 


2 


3 


1 



a Collision number, T ex p^-^r^, normalized to the value for M4, where 
po is the cluster central density and r c the core radius. Cluster parame- 
ters come from Harris (1996; updated Feb. 2003), except M4 for which 
Pq and My are computed from the distance and reddening of Richer 
et al. (1997). 

b Scaled cluster mass M^ inside r^. Following Kong et al. (2006), 
scaled values for M h were calculated from M h = l0- aA ( M v~M VtMi ) ^ 
where My is the cluster absolute visual magnitude from Harris (1996, 
updated 2003). The value My iM 4 for the cluster M4 is -7.2. 

C L X is the X-ray luminosity of individual sources in the 0.5 — 2.5 keV 
energy band 

d M4 possesses one optically identified X-ray source and likely CV with 
L x > 10 31 ergs s _1 . 

c For each cluster, the basic data for this table were extracted from: 1 
Bassa et al. 2004; 2 this work; 3 Heinke et al. 2005, but correcting for a 
distance of 4.5 kpc; 4 Grindlay et al. 2001b. 
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o 



30 



o 



4400 



Fig. 1.— ACIS-S S3 0.3-8.0 keV image of the globular cluster M71 for r Mn < 2r h . The 
small circles show the positions of the 29 X-ray sources listed in Table [1] with rM7i < r h 
and the 34 X-ray sources also listed in Table [1] with r h < r Mn < 2r h . The large circles 
are centered on the nominal center of the cluster and have radii r c (inner circle), r^ (middle 
circle), and 2r^ (outer circles). Straight lines mark the nominal boundaries of the S3 and S2 
CCDs, with most of the figure falling on S3 and the lower left hand portion on S2. 
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Fig. 2. — Number of sources, N(< tmt\), inside radius tmi\ (') vs. tmii- The top panel is 
for all S3 sources inside 2rh, and the bottom panel for those sources inside 2r^ with numbers 
of 0.3-8.0 keV counts C x > 10. The histogram shows the actual number of detected sources 
inside the corresponding radius, while the upper solid curve shows the number predicted by 
the best-fit King model. The lower curves are the model predicted background and cluster 
contributions, as indicated. The vertical lines mark the core radius, r c , and half-mass radius, 
rh- 
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N(>C X ) 
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Extragalactic 
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100 
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Fig. 3. — The M71 log iV-log C^ distribution for sources with Tmi\ < 2r^. The solid lines 
through the data represent unweighted least-squares fits to a power-law for C x (0. 3-8.0 keV) 
< 50 and 0,(0.3-8.0 keV) > 50. The functions are iV(> C x ) = 10 a / C b x , with (a, b) = (2.22, 
0.67) for C x (0.3-8.0 keV) < 50, and (a, b) = (5.10, 2.34) for 0,(0.3-8.0 keV) > 50. The 
lower lines show the estimated contribution from extragalactic sources (solid Giacconi et al. 
2001, dotted Brandt et al. 2001). 
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Fig. 4. — Color-color diagram for sources with rMn < 2r^ having 45 or more source counts 
in the energy band T = 0.3-8.0 keV, plus s08, the MSP candidate counterpart, and the 
summed spectra for Groups 1-4. Sources inside r^ are shown with large black dots while 
sources between r^ and Ir^ are shown with large grey dots. The bands are S = 0.3-0.8 keV, 
M = 0.8-2.0 keV, and H = 2.0-8.0 keV. The dot-dashed curves are the PIMMS predicted 
values for power-law spectra with indices ranging from -1 (hardest) to -6 (softest), and column 
densities n H = 1.0 x 10 20 (bottom), 1.39 x 10 21 (middle), and 1.0 x 10 22 (top) cm" 2 . The 
labeled small grey dots show the positions along these curves of the photon indices -1 to -6, 
from right (hardest) to left (softest), in increments of -1. Although we have included results 
for 52 and Group 4 on this plot, power-law models with n# fixed at the value appropriate 
for M71 did not provide acceptable fits to the X-ray spectra from these sources. 
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Fig. 5.— The 0.3—8.0 keV light curve for PSR J1953+1846A, folded at the reported radio 
binary period of 0.1768 d = 4.2431 h = 15.2751 ks. The bin size is 1/5 the binary period. 
The single period probability of chance occurrence of the corresponding value of x 2 is 2.2%. 
Also shown is the phase spanned by the radio eclipse. 
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Fig. 6. — Light curves (0.3-8.0 keV rate in c/ks vs. time in ks) for the 14 Chandra detected 
X-ray sources in Tables CQ-TS] with 0.3-8.0 keV counts C x > 80, plus the X-ray source, s08, 
coincident with the MSP in M71. Both the source and the local background light curves are 
shown. 
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Fig. 7. — For four clusters, number of background-subtracted sources inside r^ with 0.5 — 
2.5 keV X-ray luminosities, L x3 o, in the range 1 — 10, vs. scaled values for the mass, M^, 
inside r^. Following Kong et al. (2006), scaled values for Mh were calculated from M^ = 
Iq~oa(m v -m Vi m4.) ^ w \± eTe My is the cluster absolute visual magnitude from Harris (1996, 
updated 2003). The value My,M4 for the cluster M4 is -7.2. The result of a weighted linear 
fit to the data (JVj, = aMh), including M71, is also plotted, with slope a = 11.1. 



